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ABSTRACT 

We demonstrate the stabilization of an isolated unstable periodic orbit in a liquid bridge con- 
vection experiment. A model independent, nonlinear control algorithm uses temperature measure- 
ments near the liquid interface to compute control perturbations which are applied by a thermo- 
electric element. The algorithm employs a time series reconstruction of a nonlinear control surface 
in a high dimensional phase space to alter the system dynamics. 

INTRODUCTION 

Hydrodynamic instabilities arise frequently in materials processing and can significantly im- 
pact the quality of the resulting products; these instabilities can be present even in a microgravity 
environment where buoyancy effects are suppressed. In float-zone refinement, for example, the 
appearance of time-dependent thermocapillary convective flow induces undesired variation in the 
chemical composition of crystals processed either on Earth or in microgravity [1]. Recent experi- 
mental and theoretical studies have explored hydrodynamic instabilities in model systems such as 
liquid bridges [2-5]; these studies are an important first step toward devising strategies for improving 
process outcomes. 

There has been much progress in the control of dynamics in low-dimensional nonlinear systems 
[6]. These efforts, popularly known as “control of chaos”, have succeeded in stabilizing various 
physical systems including lasers [7], chemical reactions [8], and cardiac tissue [9]. Unfortunately, 
these methods fail when the goal dynamics for the controlled system is distant in phase space from 
the autonomous system behavior. Recently, however, an algorithm that targets goal dynamics 
anywhere in the phase space was proposed [10]. The algorithm requires no knowledge of the 
underlying, possibly nonlinear, governing equations; instead, time series of the system response to 
random perturbations is used to construct a control law in a time-delay space. 


EXPERIMENTS AND CONTROL ALGORITHM 

Convection Cell A liquid bridge is formed by a drop of silicone oil (tetradecamethylhexas- 
iloxane with > 95% purity [11]) confined between two coaxial cylindrical stainless steel rods (radius 
r = 0.3 cm) that are separated vertically by a gap l = 0.3 cm [Fig. 1(a)]. The rods differ in 
temperature by AT > 0 with the upper rod warmer than the lower; the mean temperature of the 
bottom rod is 15.0 °C and AT is computer-controlled to a precision of ± 0.05 °C. The volume of 
liquid in the drop is 0.065 cm 3 . The dimensionless number that characterizes the surface tension 
driving is the Marangoni number (M = a T ATdjpuK , with the liquid surface tension a, density p, 
kinematic viscosity thermal diffusivity k, and a T = |^~|). 
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The flow is measured and controlled by, respectively, a single temperature sensor and single 
feedback element. The temperature sensor is a 0.03 cm diameter thermistor that is placed approx- 
imately // 2 above the lower rod and 0.03 cm from the surface of the drop. The sensor yields a 
measurement that is determined primarily by the local temperature of the nearby liquid interface. 
The feedback element is a 0.1 cm x 0.2 cm x 0.3 cm thermoelectric device which is placed at the 
same relative location as the temperature sensor on the opposite side of the liquid bridge [Fig 1 (a)]. 
The temperature of the feedback element can be increased or decreased by applying a voltage of the 
proper polarity; these temperature changes, in turn, impose a localized temperature perturbation 
at the surface of the drop, thus altering the thermocapillary forces that drive the flow. 

Control Algorithm A representation of the system state is obtained by selected sampling of the 
temperature sensor output and the applied perturbations. Samples are recorded at local maxima of 
the sensor output; differences between adjacent maxima are taken to filter out the influence of slow 
temperature drifts. The resulting sequences of temperature differences x f and applied perturbations 
Pi (which are held constant between samples i — 1 and t) are used to construct a simplified map 
description of the dynamics through the use of difference equations. 

For an m-dimensional system at time step i, the next applied perturbation p l+ i is determined 
from the history (subscript h) and target (subscript £) dynamics: 

Phi = C [x* (i) , p h (0 , x* (i + 1) , p f (*)] 

where C is a single-valued mapping into a time-delay space with x A (t) = (x iy ..., x t ’_ m +i), p/,(i) s 
(pi — di •••? P»— m-d+2)) x$(i) = • ••, £»*+2m+d-i)t and Pt(i) = (p*+m+i i Pi+2m- 1) * This formu- 

lation assumes that measurements are recorded from a single sensor; therefore, a given system state 
can be determined by sequence of m measurements and m perturbations. Only one control param- 
eter is available to change the system dynamics; thus, a minimum of m successive perturbations 
are required to reach the target dynamics, which, for stabilization of an unstable periodic orbit, is 
given by x t (t + 1) — 0 and p*(i) = 0. Additionally, d is a delay time which includes the time to 
compute pi+i and the time before an applied perturbation affects the measurements. 

If C is linear, 4m unknown coefficients are required to determine p*+i. For the liquid bridge, C is 
nonlinear and is approximated by a tangential (linear) interpolation in the delay space constructed 
from time series. Typically, one thousand random perturbations are first applied to the liquid 
bridge. The applied perturbations and resulting measurements from this “identification stage” are 
stored as a reference set. When control is activated at time step i } the reference set is searched for 
8m readings (a minimum of 4m are required) that are nearest in the delay space to the history and 
target dynamics [x ft (i), p/,(i), 0, 0]. The best linear least squares fit to the 4m coefficients of C is 
obtained by singular value decomposition of the (overdetermined) 8m system of selected reference 
readings and the next applied perturbation for each reading. With the determination of C, the 
perturbation p l+1 is computed and applied. The entire computation procedure takes approximately 
0.1 sec on a 120 MHz Pentium PC; this forces a one iteration delay (2.3 s) before the perturbation 
is applied. 


RESULTS 

The control scheme is applied for M — 17750 where the unperturbed system exhibits a two 
frequency state. For M > ~14000, the flow becomes time-dependent with a single fundamental 
frequency [Fig. 1(b)]; for M « 16500, a second frequency appears and is locked to the first. When 
the control algorithm is applied, the second oscillation is rapidly suppressed [Fig. 2(a)]. Initially, 
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the applied perturbations are large; however, once the system approaches the target dynamics, only 
small perturbations are required to keep the system under control [Fig. 2(b)]. After the control is 
turned off, the system rapidly returns to the two-frequency dynamics of the unperturbed state. 

Effective control is obtained when C is computed for several different values of m and d ; howev- 
er, the fastest convergence is achieved for m = 4 and d = 2. This result implies the dimensionality 
of the unperturbed dynamics is m = 4. Two dimensions are required to describe the second fre- 
quency present in the unperturbed system; the other two dimensions effectively describe the decay 
of stable modes in the liquid bridge and the thermal relaxation of the feedback element. The delay 
d = 2 is approximately equal to the sum of the calculation delay and the time for waves with 
azimuthal wavenumber 2 and period of 2.3 s to propagate from the feedback element location to 
the sensor location. 

Figure 3 demonstrates that our control method is effective for stabilizing states which lie far 
from the unperturbed dynamics in the phase space. Previous methods for controlling periodic 
orbits relied on the unperturbed dynamics to carry the system close to the target dynamics before 
the application of control [6]. The target dynamics for our control method can lie in any region of 
phase space that can be accessed by the random perturbations applied during the “identification” 
period when the reference set is constructed. The control scheme fails for M > ~19000 because 
the dynamics become highly nonlinear, and the one thousand points in our reference set become 
insufficient for good interpolation. More sophisticated methods developed for nonlinear time series 
predictions neural networks) may help extend the parameter range for control [12]. 


CONCLUSIONS 

Our experiments demonstrate that a single local measurement and feedback perturbation are 
sufficient to control liquid bridge convection in parameter regimes where the dynamics are suffi- 
ciently low-dimensional. For other flow states, however, the spatial structure cannot be ignored. 
In particular, we have attempted to stabilize unstable time-independent states using the present 
experiment. Oscillations can be suppressed at the sensor location; nevertheless, infrared imaging 
reveals the presence of standing waves with antinodes between the feedback element and the sen- 
sor. In this case, multiple spatially distributed measurements and perturbations will be required 
for control; we are presently modifying our control algorithm and experiment to include two sensors 
and two feedback elements. In this way, liquid bridge convection serves as an ideal proving ground 
for developing methods for controlling spatially extended nonlinear systems. 
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Fig. 1. (a) Sketch of liquid bridge convection experiment, (b) Infrared image of the brightness tem- 
perature (darker shading for colder temperatures) for time-periodic liquid bridge convection. The 
“barber-pole” structure of the temperature corresponds to two waves that propagate azimuthally 
(right to left in the figure). 




Fig. 2. The application of control is illustrated for discretized dynamics by time series of tempera- 
ture differences ar t (top) and applied perturbations pi (bottom). The control is applied at time step 
i = 300 and the second oscillation in the two frequency state is rapidly suppressed. Releasing the 
control at time step i = 800 caused the system to rapidly return to the unperturbed dynamics. 
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Fig. 3. Second return map constructed from experimental time series illustrating both the torus 
dynamics (control “off”) and the stabilized periodic orbit (control “on”). 


492 




